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ABSTRACT 



Surface mixed layer properties off Northern California 
(NOCAL) were analyzed statistically and numerically. The 
observations were acquired on three cruises as part of the 
Pilot Ocean Prediction Study of the California Current 
eddies centered ca. 37 to 39N, 125 to 127W during March and 
August 1982. Mixed layer depth, averaging 33+14m, had a 
horizontal correlation scale of no more than 35km, which has 
significance for relating thermal structure information from 
individual temperature profiles to that of Fleet Numerical 
Oceanography Center's (FNOC) analyses based on a grid length 
of approximately 300km. Simulations and sensitivity tests 
were made using the Garwood bulk mixed layer model and the 
Mellor Level-2.5 diffusion model with the initial and 
boundary conditions acquired at sea and from FNOC. Upper 
ocean thermal structure analyses and forecasts were also 
obtained from the Navy's TOPS/TOPS-EOTS diffusion model, 
which has since become operational at FNOC. Comparisons 
of observations, analyses, and model solutions reveal 
consistent cooling and deepening by the former two models 
and excessive warming by the latter model. These 
significant differences are believed to be related to 
model resolution, model sensitivity, oceanic and atmospheric 
data quality, and spatial variability. 
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I. INTRODUCTION 



A. BACKGROUND AND OBJECTIVES 

Over the past few years, the ocean research community 
and the Navy have become increasingly aware of the need and 
measures for ocean thermal structure modeling. Not only has 
the ocean response to atmospheric forcing become of impor- 
tance because of its feedback to the atmosphere, but also 
because of the importance of ocean thermal structure in 
acoustic propagation. For these and other reasons, a grow- 
ing study has centered on the upper ocean- "that layer 
bounded by the surface which exhibits strong seasonal varia- 
tions... most notable in the local temperature profile of the 
water column." [Grabowski, et.al., 1982] The nearsurface 
structure called the mixed- or mixing layer has received 
significant emphasis. The Glossary of Meteorology defines 
this structure "as a surface layer of virtually isothermal 
water due to wind induced turbulent motion and/or free 
convection which frequently exists above the thermocline" . 

To understand the form and dynamics of the mixed layer and 
the upper ocean thermocline in general, numerous models have 
come into use. 

The objective of this research is first of all to obtain 
a description of the characteristics of the mixed layer in a 
'test block' of the ocean (ca. 37 to 39N; 125 to 127W) 
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surveyed during the first and second OPTOMA cruises (Figs. 1, 2) 
occurring during the second week of March and the first two 
weeks in August 1982, respectively. (OPTOMA- Ocean Prediction 
Through Observations, Modeling, and Analysis, -is a joint 
Harvard/NPS project ’intended to acquire field data to 
characterize synoptic scale eddies over a domain in the 
California Current off Northern California, and to 'set up' an 
eddy-resolving, statistical/dynamical, limited-domain, open 
boundary numerical ocean prediction model’.) The next 
objective was to evaluate the thermal structure analysis of 
the Navy’s TOPS model over the second cruise period in com- 
parison with the Mellor level-2.5 model and the Garwood 
(1976, 1977) bulk mixed layer model. Finally, the results of 
the three basically single point models were compared with the 
observations. (Note: TOPS is often considered a quasi- 

three-dimensional model, because it includes advection by 
climatological geostrophic currents and wind drift.) 

The idea that the vertical distribution of temperature, 
salinity, and current in the upper ocean layers generally 
appears to be governed largely by the vertical heat, momentum, 
and salt fluxes imposed by local air-sea transfer [cf. 

Camp and Elsberry, 19783 has led to one -dimensional time 
dependent models which assume horizontal homogeniety. 

[Klein, 1980] The reliance upon this assumption varies 
among models as well as regions of applicability. 
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Figure 1. California Current study domain for Mar. 

1982 OPTOMA cruise: (X)’s denote CCSI 

XBT stations 
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Figure 2. 



California Current study domain for Aug. 
1982 OPTOMA cruises: (X)'s denote CCSII 
and CCSIII XBT stations 



B. LITExRATURE REVIEW 



Techniques for modeling vary significantly in the degree 
to which they attempt to represent the upper ocean features 
and dynamics. The existing models may be categorized as 
belonging to one of two broqd classes: 

a. Bulk Mixed-Layer Models 

b. Diffusion Models 

Each class of models will be described with some reference 
to its components. No attempt, however, will be made to cover 
all of the models in great detail or to describe each of the 
existing models or model modifications. The only exceptions 
to this will be the Garwood (1977) and TOPS [Mellor and 
Durbin, 1975] models which are of direct interest in this 
study. The structure of the Mellor Level-2.5 model will not 
be specifically discussed, because it is very similar to the 
TOPS (Mellor Level-2) model. 

1 . Bulk Models 

The first category of models, the bulk models, are 
often referred to as slab models. The reason for this is 
that Pollard and Millard (1970) and Halpern (1974) showed, 
with the occurrence of strong wind-forced vertical mixing, 
the mixed layer responds at the inertial frequency as a 
rigid slab. This description may be misleading, for it 
implies that T, S, and V are uniform throughout the layer. 
(Glossary of Terms) Such, however, is not the case. Instead, 
the models are based on the assumption that values of the 
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aforementioned parameters are 'quasi-uniform' within a 
'well-mixed layer'. [Klein, 1980] Provided the assumption 
holds (as it may away from fronts and eddies), the values 
can be 'lumped' into integrals. [Niiler and Kraus, 1977; 
Zilitinkevitch et.al, 1979] That is, the 'bulk' models 
predict integral (or average mixed layer) values for all 
variables. (Garwood, personal communication) -hence , the 
perhaps 'better' term -'integral model' is frequently used. 

This type of modeling requires an a priori assumption 
that a mixed layer exists. [Marchuk, et.al., 1977] To 
describe the layer requires only knowledge of the fluxes 
at the surface and the base of the mixing layer (i.e., 
surface and entrainment fluxes). The mathematics of the 
modeling is, therefore, greatly simplified from partial 
differential equations to ordinary differential equations. 

Basically, mixed layer modeling is founded on two 
fundamental hypotheses. The first is that the mixed layer is 
formed as a result of the upper ocean response to forced and 
free convection. Atmospheric forcing in the form of applied 
wind stress converts mean kinetic and potential energy to 
turbulent kinetic energy within the upper ocean. This 
turbulent motion, in turn, entrains below -layer water into the 
layer, thereby changing the layer's characteristics. Strong 
solar heating works against the forced convection to 
increase the potential energy, and form a shallow turbulent 
boundary layer. This heating in effect alters only the near 
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surface water, for invoking the Second Law of Thermodynamics, 
water cannot be unmixed. Thus, solar heating results in a 
new mixed layer in addition to the ’old mixed layer’. 

(Fig. 3) This distinction between new and old mixed layers 
will be discussed and utilized in the following sections. 

The second hypothesis of importance, a topic of much 
controversy, invokes the use of some form of the mechanical 
energy budget as essential in the closure scheme of the more 
recent models. The disagreement, and consequently a differ- 
ence in some models, lies in the question as to whether the 
kinetic energy of the mean flow or the turbulent kinetic 
energy plays the main role in driving the downward buoyancy 
flux at the base of the turbulent boundary layer, i.e., in 
controlling the mixed layer dynamics and entrainment. 

Thompson (1976), Pollard (1977), and Price (1977) advocate 
that the mean kinetic energy is converted directly to poten- 
tial energy by an instability in the mean flow that causes 
the bulk Richardson number to be less than some constant of 
order one. Advocates of the turbulent kinetic energy 
budget, following the original argument of Kraus and Turner 
(1967), disagree that this will be the dominant mechanism, 
reasoning that turbulence-generated instabilities normally 
have the capacity to erode the interface to the extent that 
they preclude the existence of mean flow instabilities. 
[Garwood, 1979] 
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h, • mixed layer depth ( summer) 
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Figure 3. Temperature versus depth profile with old 
and new MLD's; h^=new mixed layer depth, 
h„=old mixed layer depth [Husby and 
Nelson, 1982] 
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Differences in the views by no means end there; for 
example, there is no unanimous approach to parameterization 
of the kinetic energy budget. The gross parameterization 
approach is commonly used in closure schemes. However, 
second order closure schemes have been developed, for example 
Launder (197 5). Such closure schemes, though applied in only 
one of the bulk models - Garwood (1976, 1977) -are receiving 
increasing use. Garwood's scheme will be described in a 
later section. 

2 . Summary of Bulk Models 

The structure of all simple one-dimensional bulk 
models may be described in terms of two fundamental equations 
the vertically integrated heat equation and the vertically 
integrated turbulent kinetic energy (TKE) budget. 

(1) h(3T /at) = Q - 0, 

o o h 

A B C 

h h 

(2) P = / Gdz+E -E,+6( f Fdz- i hFl__, ) 

Q oho ^ z-n 

D E F G ^1 ^2 

^ ^ ae 

- / edz - / dz. 

o o 

I J 
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where indicated terms are: 

A: Mixed layer depth (MLD) times the local time rate of 

change of mixed layer temperature (MLT). 

B: Net heat flux through the sea surface. 

C: Net heat flux through the base of the ML. 

D: P= iehCQ^ + Q^) 

Vertically integrated generation or dissipation of 
TKE by buoyancy forces. 

E: Rate of TKE production due to transfer from mean 

current shear. 

F : E^ is the downward TKE flux at the top of the ML 

(i.e., from the atmosphere to the ocean and 
produced by wave breaking at the sea surface). 

G: E^ is the downward turbulent energy flux at the base 

of the ML. 

and H 2 : Effect of penetrative solar radiation. 

I: Net turbulent layer dissipation over the ML. 

J: Net local time rate of change of TKE over the ML. 

The main differences between the various bulk models 
are in the definitions of the heat flux, Q, and the turbu- 
lent entrainment across the interface. Some of the main 
features of several bulk models are listed below. 

Kraus and Turner (KT) (1967) ; KT developed the first bulk 
. model capable of addressing the time evolution of a mixed 
layer/seasonal thermocline system; they were the first to 
propose a TKE budget to close the system. Kraus and Turner 
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assumed a steady state turbulent regime and neglected 
dissipation. They parameterized penetrative solar radiation 
with an exponential decay law and shear production in terms 
of the frictional velocity. 

Denman (DN) (1973) : DN rederived the KT model, following 

Kraus' advice that the rate of increase of potential energy 
by mixing should be a small constant fraction, m, of the 
rate at which turbulent kinetic energy is transferred down- 
ward by the wind stress. (He used m=.0012) Denman also 
included absorption of solar radiation below the mixed 
layer. He showed that for sufficiently long model runs, the 
eventual depth of the mixed layer depended as much upon the 
temperature gradient below the mixed layer as upon the 
generation of turbulent kinetic energy. Unlike KT, DN 
included some turbulent dissipation. 

Pollard, et.al. (PRT) (1973) : As DN and KT , PRT assumed a 

steady state regime. The rate of work done by the wind 
stress on the mean motion was equated with the rate of 
increase of the energy of the mean field. The mixed layer 
deepened rapidly to a maximum within half of an inertial 
period from the time of the instantaneous increase in the 
wind speed (which then was held constant). Forty percent 
of the incoming energy was found to directly increase the 
mean kinetic energy while the rate of change of the mean 
potential energy was three times smaller than that of the 
mean kinetic energy. This result reinforced their 
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hypothesis that the shear instability of non-steady 'wind 
drift currents' is the main mechanism in turbulence 
production. They neglected turbulent dissipation and 
penetrative solar radiation. 

Niiler (NI *75) (1975) : NI recognized that instantaneous 

sharp wind changes were not realistic, and that in a 
steady regime, the surface energy fluxes from the breaking 
waves could be an order of magnitude greater than the energy 
fluxes from the wind to the currents. Thus, NI extended the 
PRT model so that at the beginning and final stages of a 
model run, TKE generation through surface fluxes was most 
important, while at intermediate times the drift current 
velocity shear was dominant. As with the PRT, NI'75 
assumed a steady state, no penetrative solar radiation, and 
no energy flux through the base of the mixed layer. 

Turbulent dissipation was parameterized as a function of 
the frictional velocity. 

Kim ( KIM) (1976) : KIM proposed a model more complex than KT 

which included TKE storage as well as the same parameteri- 
zation for penetrative solar radiation as in DN. An 
additional feature of importance was the use of a background, 
depth-dependent turbulent dissipation rate to avoid infinite 
deepening . 

DeSzoeke and Rhjnes (DR) (1976) : They used the NI model 

with the addition of a TKE storage term (term J) to show 
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that the KT and PRT models yielded results that are 
assymptotes for a single model equation. 

Gill and Turner (GT) (1976) : They simulated oceanic 

seasonal cycles, and showed that the cyclical behavior of 
the potential energy content, heat content, and the 
surface temperature could be recovered if all or some 
portion of the potential energy released during the runs 
was dissipated. They neglected terms G, H^, H2, and J. 

Niiler (NT ’77) (1977) : This model is similar to NI '75. 

Additions include a TKE storage term (term J) and background 
dissipation identical to that employed by Kim. 

Niiler and Kraus (NK) (1977) : The NK model is a conglomerate 

of several of the models described above. For example, the 
radiation term is the same as that used in the KT model; 
terms E and F are similar to those used previously by 
Niiler with the only difference being in the coefficients; 
and the turbulent dissipation term is similar to that used 
by Resnyansky (1975) and the shear production term of the GT 
model. The NK model does not include the storage term, J. 
Unlike previous models, it includes the downward flux of TKE 
through the mixed layer base (term G). 

Garwood (GWD) (1977): (to be discusse-d in Section 3) 
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TABLE I** 

SUMMARY OF TERMS USED IN BULK MODELS 
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KT 
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+ * 




DN 


* 






PRT 
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NI '75 




•*: * 




KIM 
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DR 


- 




* 


GT 








NI '77 
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4% 4% 


* 


NK 


* 


SS 4^ 4\ 




GWD 




- * * 




" Appears explicitly 

+ Appears implicitly, i,e., in one of the other terms 
** Modified (Zilitinkevich. et . al . , 1979) 




3. 


Garwood Model 








The Garwood (1977) bulk 


mixed layer model is 





composed of a closed system of seven equations: 
The entrainment flux equation: 



(3) 
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The horizontal component of the turbulent kinetic energy 
equation : 



(4) 
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The vertical component of the turbulent kinetic energy 
equation : 



(5) 
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The mean buoyancy and complex mean momentum equations : 



(6) h = bw( -h) -bw(0 ) + Q 
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The jump conditions (relating entrainment fluxes to mean 
buoyancy, rate of deepening, and mean momentum): 



(8) -bw(-h)=AB . 

0 L 



(9) -cw(-h)=AC , 

a T 



Depending on which version of the model is used, 
one needs either to provide surface momentum and buoyancy 
fluxes, mean buoyancy and momentum below the mixed layer 
along with the radiation absorption or else the observed 
variables from which these fluxes can be computed. Closure 
is achieved with the bulk buoyancy and momentum equations 
and the vertically-integrated individual turbulent kinetic 
energy components. 

A number of significant differences exist between 
this model and other bulk mixed layer models that may not be 
obvious from the preceding summarization or the matrix of 
terms describing which terms are included in each model. 

Some of these will be made explicit in the following 
discussion of the mechanics of, and ideas beind, the Garwood 
model. 

According to Garwood (1977), the basis for the 
initial interface destabilization and the resulting entrainment 
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is the shear across the interface (base of the mixed layer). 
This shear is hypothesized by Garwood to be the result of 
local turbulent eddies (wind generated turbulent kinetic 
energy available for mixing is explicitly parameterized as 
dependent on h/L) with the mean flow only making a minor 
contribution to the achievement of this critical shear 
value. (The instability related to this criterion is often 
called a Kelvin-Helmholtz instability or a Benjamin (1963) 
Class C instability.) By using the familiar equation for 
the flux Richardson number: 

(10) Rf=bw/(uw3U/3z+^3V/3z) , and 

taking into account the jump conditions (8, 9), shear 
production can be shown to be only a fixed fraction of the 
buoyant damping in the entrainment zone. This- zone may in 
fact have a Rf> 1, and still possess enough turbulent 
kinetic energy for entrainment to proceed. Garwood’s 
reasoning, in contrast to that of Pollard, et ,al . (1973), 
is that the mean shear production is only a secondary 
energy source for mixing. This secondary source of energy 
can only be made available by entrainment which is 
initiated by another source. This other and more 
significant source is the convergence of the flux of 
turbulent energy from above by the term: 

(11) |-[w(|-+ ^)] 

3 Z p Q 
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Therefore, the critical number in the Garwood model is not 
the flux Richardson number, but instead the ratio of the 
buoyancy flux to the convergence of turbulent energy, P: 



( 12 ) 




The time, x 



e 



(13) 



a^h<w > 



- 1/2 



X 



e 



required to transport the available turbulent energy, <E> , 
to the entrainment zone is important in relation to dissipa- 
tion, and to the idea that not all available TKE contributes 
to increasing the potential energy of the system. According 
to Tennekes and Lumley (1972), the viscous dissipation of 
geophysical flows with large Reynolds numbers is proportional 
to the reciprocal of the time scale of the largest eddies. 
Garwood developed the idea that for the deeper mixed layers 
(R **1), the planetary rotation introduces a second time scale 



of importance: (1/f). Not only did the inclusion of this 

second time scale enable Garwood to more explicitly 
parameterize viscous dissipation in relation to the local 
Rossby number, but the very role of the time scale in the 
entrainment equation enabled Garwood’s model to simulate a 
wider range of conditions, i.e., diurnal and annual ranges 
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of mixed layer stability. (The PRT model can simulate a wide 
range of changes, though the deepening of the mixed layer is 
in accordance only with the mean shear production.) 

Another difference over most previous models was in 
Garwood’s method of using the turbulent kinetic energy 
equation. The model uses the component equations (4, 5) 
vice the typical use of the total equation. This facilitates 
analyzing the retreat or shallowing of the turbulent boundary 
layer which depends on only the vertical component of turbu- 
lence being inadequate to transport heat, momentum, and TKE 
to the earlier depth of mixing. (Fig. 3) Thus, using the 
equation in component form allows for a more explicit and 
satisfactory treatment of the mixing process. 

4 . Diffusion Models 

The second category of mixed layer models, the dif- 
fusion models, has received less attention from the ocean- 
ographic community than the simpler bulk models. Yet, the 
interest that exists has grown greatly with the increasingly 
powerful computers and the development of methods for solving 
nonlinear partial differential equations. [Marchuk, et.al 
1977] The diffusion models are both complex to interpret 
and to program. They require the numerdcal solution of 
P'artial differential equations of temperature, salinity, and 
momentum in time and depth for closely spaced vertical 
arrays of gridpoints. Unlike bulk models, the 'gridpoint 
models’ do not require the a priori assumption of the 
existence of a mixed layer. Mixed layer depths, are, on the 
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contrary, diagnosed from the examination of the turbulence 
flux predictions as well as the predicted temperature and 
salinity profiles. 

Mathematically, the ’mixed layer' in a diffusion 
model's vertical profile is treated in the same manner as 
is the thermocline. This allows for a more detailed analysis 
of the vertical structure (mixed layers, thermoclines , etc.), 
the mixing processes throughout the water column, and a more 
complete simulation of the vertical fluxes of temperature, 
salinity, and momentum. The latter is a basic concern in 
the diffusion models which attempt to represent the turbulent 
fluxes w ' u ' , w ' T ' , and w * S ' by algebraic or differential 
equations in terms of the gradients of mean temperature, 
salinity, and momentum. These fluxes in conjunction with 
the necessary boundary conditions for T, S and u provide a 
solution for the basic conservation equations of temperature, 
salinity, and momentum. Thus, the system of turbulence and 
thermodynamic equations are closed in order to solve for 
the vertical distributions of thermodynamic and turbulence 
variables. [Grabowski, 1980; Kondo , et.al., 1978] 

Diffusion models may be further categorized by 
complexity, i.e., a) no-equation models; b) one-equation or 
k-models ; c) two equation or k-l-models ; and d) second 
order closure models. 

A no-equation model is based on the turbulence 
transport being parameterized by a buoyancy influenced eddy 
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viscosity coefficient. The one-equation or k-models are 
similar to the no-equation models though they include a 
relationship between eddy diffusion and the local turbulence. 
The two-equation or k-l-models include a transport differen- 
tial equation. The second order closure models, first 
proposed by Chou [1945] , are based on transport equations 
for all of the Reynolds tresses derived from the Navier 
Stokes equations. They also include turbulent heat and 
salinity fluxes and the variance derived from the latter in 
conjunction with the heat conservation equation. 

A detailed profile is sought to represent the 
oceanic vertical structure as completely as possible. For 
example, Halpern (1976) among others, noted the existence of 
shear in the top of the water column, a characteristic 
assumed not to exist in bulk models. Subsequently, the 
necessity of using numerical models to determine a more 
realistic distribution of variables such as the velocity, 
temperature, turbulence intensity, dissipation, and the like 
was recognized [Phillips, 1977]. 

The first ocean diffusion model was implemented by 
Ekman in 1905. It was a type A model, that is, it was a 
no-equation model that neglected buoyancy and salinity. A 
notable result from this model was the famous Ekman spiral. 
Because of its neglect of density stratification, this model 
will not be discussed further. The major difference between 
the diffusion models, beyond the number of equations 
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employed, lies in whether they are based on observations, 
dimensional analysis, and simple stability arguments or on 
some form of the TKE energy equation. Of the five models to 
be briefly discussed below, the first two are of the former 
type, while the other three are of the latter type. 

5 . Summary of Diffusion Models 

Munk and Anderson (MA) : 1949: The HA model is a 

no-equation model wherein the eddy diffusion coefficients 
for heat and momentum are functions of the gradient 
Richardson number. No salinity is included in this model. 

Kondo, Sasano, and Ishii (KSI) : 1979: The main fea- 

ture of the SKI model is that constant flux formulations are 
applied to regions of varying fluxes. Formulation of 
diffusion coefficients based on studies of the atmospheric 
boundary layer are cast in Monin Obukhov formalism and 
applied to the upper oceans. Eddy diffusion coefficients 
for salt are assumed to be the same as those for heat. 

Vager and Zilitinkevich (VZ) ; 1968: This model is 

a one-equation or k-model based on Prandtl’s extension of 
the mixing length theory for homogeneous flow. Salinity is 
ignored . 

Mellor and Yamada (MY) : 1974: Mellor and Yamada 

formulated a sequence of models with varying degrees of 
complexity using second order closure. This discussion will 
deal with the level-2 model, because the TOPS/TOPS-EOTS 
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model (to be discussed in the next section) is based on the 
level-2 MY turbulence closure scheme, and because it is simpler 
than the level-2.5 model. (.The level-2.5 MY model was also 
used though its mechanics are quite similar to those of the 
level-2 MY model. The main difference is the inclusion of 
additional complexity as, for instance, the retention of the 
TKE storage term in the level-2.5 MY model.) (To avoid 
repetition, TOPS/TOPS-EOTS will not be mentioned in this 
section except to say that the ’Clancy version’ employed 
includes salinity, solar radiation, and horizontal advection 
by instantaneous wind drift and cliraatologically averaged 
geostrophic currents.) The level-2 MY model neglects 
horizontal advection and diffusion terms. The TKE budget is 
consequently a balance between the production and dissipation 
of TKE and the conversion to potential energy as seen in 
Table II. A. similar balance exists between the turbulent 
heat flux, the temperature variance, and the individual 
stress components. Unlike the other two models described in 
Table II, the MY level-2 model neglects the turbulent transport 
term. 

Kochergin, Klimok , and Sukhoruks (KKS): 1976; The 

KKS model is a two equation or a k-l-model. Its TKE budget 
parameterizes the time rates of change of TKE and the 
turbulence dissipation rate. 
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TABLE II 



SUMMARY OF TERMS USED IN DIFFUSION MODELS 
(Modified: Grabowski, et.al. 1980) 



TKE Prod/ Shear 
dissip Produc- 

rate tion 



Buoyancy 
Conver- 
sion of 
PE to TKE 



Vertical 
Transport of 
TKE Work Done 
against Pres- 
sure Gradient 



Dissipation 



VZ 



0 



MY 0 



0 0 



KKS 



0 



includes term 

0: does not include term explicitly 



6. TOPS/TOPS-EOTS 

TOPS /TOPS-EOTS is the Navy's new real-time ocean 
thermal structure analysis/forecast system in use at the 
Fleet Numerical Oceanography Center (FNOC). (Operational as 
of March 1983) It was developed by the Naval Ocean 
Research and Development Activity (NORDA) as part of the Navy's 
automated environmental prediction system (AEPS). The goal 
behind TEOTS is to produce reliable, accurate, real-time 
representations of the ocean thermal structure. The motivation 
for the way the model is organized (to be discussed below) lies 
in the available sources of data and the operational nature 



41 



of the model. Given that the upper ocean is primarily 
atmospherically forced, and the assessment that the mixed 
layer depths have proved ’highly predictable with a variety 
of models', such as, the Denman (1973) or Mellor and Durbin 
(MD) (1975) models, NORDA has developed an operational model 
based largely on the use of the rather well-defined atmos- 
pheric variables to drive a model to describe the relatively 
data sparse oceans. Surface wind, solar heat fluxes, and 
precipitation fields are among those necessary to drive 
TOPS/TEOTS. 

For vertical structure, TOPS has seventeen levels 
from the surface to 500 meters, with the majority of the 
levels placed within the first 100m to enhance the vertical 
resolution where most needed. (Appendix B) The horizontal 
grid is typical of many operational atmospheric models, 
i.e., 381 km at SON. The input fields on this grid are 
smoothed spatially by the TOPS modified EOTS . Spatial 
averaging of the fields also smooths the temporal varability 
of the fields. Of course, horizontal features smaller than 
the spacing cannot be resolved so that the capabilities are 
limited in enclosed seas, in the vicinity of fronts and 
eddies, and near boundary currents. [Clancy and Poliak, 1982] 
This last point is important to keep in mind for the 
simulation comparisons in Chapter IV. 
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As its name implies, TOPS/TOPS-EOTS is made up of 
two components: TOPS and TOPS-EOTS. They are coupled in a 

cyclical fashion. 

TOPS-EOTS is an improved, updated version of the 
Navy's conventional objective analysis component, EOTS 
(Expanded Ocean Thermal Structure). It supplies TOPS with 
updated (real-time) initial conditions for a twenty-four 
hour forecast. Its improvement over conventional EOTS lies 
in the fact that the TOPS forecast is fed back into TOPS-EOTS 
as a first guess field for the following day's analysis. This 
not only reduces the noise resulting from the introduction of 
data into the analysis, and makes the analysis dynamically 
and thermodynamically consistent with the atmospheric 
forcing, but it also allows the model to be run for an 
extended period of time within data-sparse regions prior to 
being updated, provided atmospheric forcing data are available . 

TOPS (Thermodynamic Ocean Prediction System) , the 
forecast component, is a thermodynamic model for vertical 
structure which is forced by surface fluxes supplied by 
NOGAPS (Naval Operational Global Atmospheric Prediction 
System). It uses the MY level-2 turbulence scheme previously 
described with the added effects of advection by 
instantaneous wind drift and climatologically averaged 
geostrophic currents. 
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The prognostic equations used by the model are the 
conservation equations of temperature, salinity and momentum 
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where : 

*D' is a damping coefficient representing the drag force by 
radiational stress at the base of the mixed layer associated 
with the propagation of internal waves downward and away 
from the wind-forced region. ’A’ is the background 
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horizontal eddy diffusion associated with the intermittent 
breaking of internal waves. [Clancy, et.al., 1982] 

There are no horizontal pressure gradients in the 
last two equations so that the u and v represent wind drift 
current. The geostrophic advection is provided to the u 

“a 

and v^ terms in the T and S equations. 

Within the MY level-2 closure scheme, the vertical 
eddy fluxes of T, S, and momentum are parameterized 
(Mellor and Durbin (1975); Clancy and Poliak (1983)) as 
follows : 

(17) (T-r = -lqs„ II H -K„ 

(18) -K„|| 
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where K„ and K„ are eddy diffusion coefficients; 5. is a 
M H 

turbulence length scale; q=sqrt ( ( 2 ) *TKE) ; and are 
functions of the bulk Richardson number: 
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There is an implied cutoff Richardson number such 
that one can determine the three empirical constants. The 
TKE budget used to obtain q is: 
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and the equation for I is: 
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The solar radiation flux is parameterized by using an 
extinction profile for the most common water type. (For 
the level-2.5 model, the solar radiation flux is applied 
at the surface.) 

As mentioned previously, boundary conditions are 
mainly provided by the FNOC atmospheric models. XBT profiles 
are included in the analysis when they are available. (The 
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ocean data used in the TOPS runs will be briefly mentioned 
in later sections when the model results are discussed.) 

7 . Verification Testing of Mixed Layer Models 

Numerous tests have been completed or are in the 
planning stages for verification of the forecasting 
abilities of the Garwood, TOPS, and Level-2 Mellor or Mellor- 
Durbin (MD) models as well as for the other models pre- 
viously discussed. A review of any mixed layer model 
verification testing would not be complete without at least 
mention of MILE, the Mixed Layer Experiment, which consisted 
of an intensive examination of the upper ocean in xhe environs 
of Ocean Station P during a 20-day period in autumn 1977. 
Excepting that brief discussion, the mainstay of the review 
of verification testing will encompass the three models 
directly involved in this study. (Though the MD model 
is essentially the core of the TOPS model, each will be 
reviewed separately, for TOPS, a first generation operational 
model, has its own inherent verification problems.) 

Mixed Layer Experiment, 1977 : The Mixed Layer Experiment 

(MILE), 1977, is a noteworthy example wherein data were 
acquired, and compared initially against sim.ulations produced 
using the NI’77 model and subsequently other model simulations. 
Davis, et.al. (1981) used velocity profiles, measurements not 
usually available in ML verification studies in an attempt 
to not only obtain a firm understanding of the dynamics and 
thermodynamics involved, but in so doing, to also confirm 
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the applicability of the bulk approach to the study domain. 
They failed to acquire such a confirmation in all cases 
except those involving strong storms wherein the upper 
oceanic shear became localized a short distance below the 
isothermal layer. However, by accounting for the vertical 
advection and horizontal pressure gradient within the domain, 
Davis, et.al. (1981) were able to indicate a ’reasonable 
balance’ in both the heat and momentum budgets, respectively. 
Two major difficulties encountered were a) the need to 
acquire accurate data representing changes below the mixed 
layer so as to better understand the vertical distribution 
of heat and momentum within the mixed layer, and b) the need 
to acquire better expertise in being able to separate effects 
of internal waves, energetic inertial motions and internal 
tides from responses of the upper ocean to storms. [Davis, 
et.al., 1981] Despite the numerous weaknesses in the study 
described by Davis, et.al. (1981), simulations encouragingly 
seem to agree with the observations. Their ’most favorable’ 
result (Fig. 4) was for m^ = 0.39 and m 2 = 0.48, where m^ 
and m 2 are tuned coefficients determining ’limits of 
excursion’ and ’rate of response’, respectively. 

a. Review of Verification Testing of Three Models 
The Level-2 Mellor model and similar models have 
been evaluated in numerous and varied simulation experiments. 
Mellor and Durbin (1975) obtained favorable results from 
tests using five weeks of data acquired at Ocean Station Papa 
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Figure 4. MILE observations versus simulations: 

NI *77, Dashed curves: observed hour 

average T(5m) and depth where T is O.IC 
colder than T(5m); solid curves; 
modeled values of mixed layer temperature 
and depth [David, et.al., 1981] 
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(SON, 145W) , as well as that data acquired nearby during 
MILE. The model's ability to handle the response of the 
upper ocean to an atmospheric low pressure system, i.e., 
the observed versus the predicted SST during and after 
frontal passage, was explored by Price et.al (1978) for the 
Gulf of Mexico region. The model's ability to amply predict 
the ocean’s response to stronger atmospheric forcing was 
verified in the case of the passage of hurricane Eloise by 
the fortunate existence of a NOAA data buoy, EBIO, located 
on the track of the hurricane. Klein (1980) used the MD 
model to simulate the variability of the mixed layer depths 
in the Mediterranean Sea. The model favorably simulated 
(Fig. 5) [Clancy, et.al., 1981] the diurnal mixed layer 
response observed during BOMEX. (The Barbados Oceanographic 
and Meteorological Experiment.) 

The Garwood model has also undergone, and is still 
undergoing, verification testing, though by a more limited 
populus. One of the most notable simulations extended over 
17 years of Ocean Station Papa data. Not only was this the 
first simulation of its kind in length, and in the handling 
of interannual variability in mixed layer response, but it 
was also the first simulation conducted for the period of 
the spring transition. (Comparison of the hindcasts with 
the observations lends encouragement to the ability of one- 
dimensional (bulk) models to account for a large part of the 
variance on time scales from the synoptic to greater than a 
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the BOMEX Experiment [Clancy, 



year. [Garwood and Adamec , 1982] Additional verification 
testing includes various studies of model performance along 
the shipping track between San Francisco and Hawaii [Steiner, 
1981], and during POLEX, a component of the North Pacific 
Experiment. [Shook, 1980] Data assimilation tests were 
performed by Warrenfeltz (1980) and by Larsen (1981). 

Because the MD model is at the core of TOPS, one 
might say that the TOPS model has undergone extensive testing 
and evaluation. Such, however, is not the case in that 
TOPS is a first generation operational model, and as such is 
subject to inherent problems not experienced by other one- 
dimensional or quasi-three-dimensional m.odels. The TOPS 
verification is, for the most part, a test of the quality of 
the FNOC operational data base, the source of not only the 
initial and upper boundary conditions, but also the verifi- 
cation fields. These verification fields have necessarily 
been extracted from the objective analysis or TOPS modified 
objective analysis scheme, for observations have rarely been 
available near the model’s widely spaced gridpoints. The 
distribution of observations is neither uniform nor 
constant. The model is, therefore, subject to not only a wide 
range of testing conditions, but also to a greatly fluctua- 
ting source of observations. When considering the advective 
version of TOPS, one must also consider possible 
three-dimensional processes. 
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Thus, meaningful quantitative forecast verification 
seems reasonable only in regions containing an adequate 
supply of data. The verification fields must be obtained from 
the same scheme used to acquire initial conditions, because 
different schemes applied to the same data set are likely to 
produce varying verification fields. [Barnett, et.al., 1980] 

The present FNOC daily global intake of data consists of 
approximately 200 XBT ’ s , 2,000 'bucket' temperatures, and 
approximately 20,000 satellite SST's from TIROS-N, This 
apparent lack of sufficient subsurface definition has 
resulted in much of the verification effort centering on the 
SST. 

TOPS is an operational model, one component of a 
resource-limited naval facility model network. As such, it 
neither has number one priority in the FNOC objectives list, 
nor can it avoid being impacted by the ever evolving compu- 
tational environment. The rate of personnel change, in 
combination with the shortage of oceanographers on the staff, 
has resulted in most of the verification testing being left 
as a milestone to be achieved in the future. 

One of the major changes which has recently 
occurred at FNOC, and which will alter the workings of TOPS, 
is the replacement of the old atmospheric primitive equation 
(PE) model with the Naval Operational Global Atmospheric 
Prediction System (NOGAPS). This change took place on 
3 August 1982, midway during the model simulation runs of 
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interest in this thesis. The change was made to improve the 
higher priority atmospheric products supplied to the fleet. 

A study of the effects of NOGAPS on TOPS results is in the 
plans for the near future. It is hoped that NOGAPS will 
provide more accurate heat fluxes for ocean prediction 
analysis . 

The analysis system has evolved from OTS (Ocean 
Thermal Structure) analysis system to EOTS (Expanded Ocean 
Thermal Structure) analysis system and finally to TEOTS or 
TOPS-EOTS (TOPS modified EOTS). The last change in the 
evolutionary sequence was implemented in full on 22 March 1983- 
that is , TEOTS replaced EOTS as a bottom boundary condition 
for the atmospheric models. Previously, EOTS was the bottom 
boundary condition of the atmospheric models so that there 
was only a one-way influence between the atmospheric and 
oceanic models, i.e., the atmospheric model drove the 
oceanic model. With the implementation of TEOTS, the feed- 
back loop is complete. Thus, the bottom boundary condition 
of the atmospheric model is more influenced by the physics 
of the situation. 

Another change to be more fully discussed later 
was a more resilient return of the model to climatology. This 
change took place on 20 September 1982. [Pollak, personal 
communication] Though not yet documented, the change 
supposedly results in the model forecast of SST remaining 
within approximately 3C of climatology with any deviations 
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being suppressed within approximately a two-week period. (The 
scientific basis for this criterion is unknown.) (The 
climatology receives a weight in the process such that sharp 
deviations from climatology are reduced by an increase of the 
weighting applied to the climatological factor.) 

Most of the published verification reports describe 
TOPS model performance prior to these aforementioned changes. 
For example, VJarn Varnas, et.al. (1982) performed a 60-day 
simulation in the TRANSPAC region of the central North Pacific 
for November and December 1976. This was done using TOPS in 
conjunction with OTS. Comparison of SST results from both 
the advective and nonadvective versions of the model revealed 
consistent skill over both climatology and persistence. 

(Fig. 6) (As will be discussed later, results also revealed 
a heat flux bias confirmed previously by comparisons of the 
heat fluxes with those calculated by Nate Clark.) Clancy, 
et.al., (1981) performed a short time scale synoptic veri- 
fication of TOPS/EOTS (nonadvective version on Cyber 175) 
using both pattern-of-change correlation techniques 
[Dobryshman, 1972], and the apparent forecast error 
technique. The former techniques demonstrated TOPS* ability 
to routinely produce real-time forecasts of large scale 
changes of SST. The latter technique indicated that, with 
adequate data coverage and limited domains of strong 
atmospheric forcing, the model could achieve a skill twice 
that of persistence. 
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The latest verifications of the coupled TOPS 
include TOPS/TEOTS vice TOPS linked with EOTS. TEOTS has 
proved to be less noisy and more consistent with the 
predictions of TOPS and atmospheric forcing; this result 
is especially true in relation to the spring transition. 
Composites of forecast verification statistics for May 1981 
indicate the skill of TOPS in forecasting both MLD and SST 
change out to 72 hours. RMS errors indicate TOPS betters 
persistence out to (but not including) 72 hours for MLD, 
but only out to 24 hours for SST. The latter deficiency 
is believed to be the result of the possible heat bias of 
the FNOC primitive equation model. It is yet to be 
determined whether NOGAPS can resolve this problem. 

C. STUDY AREA AND ITS CLIMATOLOGY 

The domain investigated is (Figs. 1, 2) located in the 
California Current System, just south of the Mendocino 
Escarpment (Fig. 7) in the region from approximately 37.5 to 
39. 6N and 125 to 127W. This system is composed of four 
currents: the California Current, the Davidson Current, the 

California Undercurrent, and the Southern California Current. 
Of main influence in the designated domain is the California 
Current, the eastern boundary current of the large anti- 
cyclonic gyre centered near the Hawaiian Islands. The 
current, an extension of the Westwind Drift (Sverdrup, 
et.al., 1942), occurs between the North Pacific atmospheric 
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Figure 7. Study domain in relation to bathymetry off 
Northern and Central California (Contour 
Interval: 200m) 
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high system and the semipermanent thermal low positioned 
over central California. Geostrophic and surface Ekman 
components of the flow are consistent with the temperature 
distribution, and the atmospheric high pressure system with 
its associated wind stress. An anticyclonic eddy is the 
major feature depicted by the dynamic topography computed 
for the CCSI domain. The dynamic topography for CCSII 
clearly depicts the anticyclonic and cyclonic eddies toward 
the northeast and southwest regions of the study domain which 
dominated the region during that time period and into CCSIII. 
[Reid, et.al, 1958] (Figs. 8-10) Nelson (1977) compiled a 
climatological description of the surface wind stress and 
wind stress curl over the California Current System, by 
computing the wind stress using the bulk formula in one 
degree squares from marine weather observations, primarily 
ship reports, archived by the National Climatic Center. 

(Figs. 11-13) 

Nelson encountered many of the typical problems that 
arise when dealing with multi-source data as, for example, 
uneven and dissimilar sampling in space and time. Only 
approximately 12% of the observations were actually measured. 
The resulting values provide a climatological view of the 
area. The domain of interest falls v;ithin a region of 
monthly mean maximum wind stress values for August, 
(approximately 1 dyne /cm* *2) The upper ocean is undoubtedly 
influenced by such forcing. Past evidence has indicated 
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CalCOFI map of mean August geostrophic flow 
at the surface for years 1950-1965 ( Dyn . m: 
interval: 0.04m) 
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50 km 



Figure 9 . Dynamic topogra,phy of study domain during 
CCSI (Dyn. cm, interval 1.0, surface 
relative to 450m) 
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50 km 



Figure 10. Dynamic topography of study domain during 
CCSII (Dyn, cm, interval 1.0, surface 
relative to 450m) 
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NCC - TCF-n - OC OECREE SUrVWIZRTION 




Figure 11. Distribution of observations per 1-degree 
square; the contour interval is 2,500 
observations; Values greater than 5,000 
are shaded [Nelson, 1977] 
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Figure 12. 



Long term mean August 
vectors: 1857-1972. 



wind stress 
[Nelson, 1977] 
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Figure 13. 



Long term mean August wind stress curl 
field: 1857-1972 [Nelson, 1977] 
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that the responsiveness of the coastal regime to such 
forcing. It is speculated that elements of the coastal 
response may be advected into offshore regions (and vice 
versa) by the predominance of jets within the region. (The 
small negative value for wind stress curl indicates slight 
open ocean downwelling within the study domain itself, 
whereas the coastal regions indicate strong upwelling 
tendencies.) Another point of speculation is the direct 
forcing of the subsurface eddy field by the wind stress curl 
mechanism [Mooers, personal communication, 1983]. 

Thermodynamical atmospheric forcing fields are also 
important to mixed layer evolution. Nelson and Husby (1981) 
compiled monthly mean heat flux fields over the California 
Current region. As with the mean wind fields, the monthly 
mean radiative and turbulent heat fluxes were computed from 
archived surface marine weather reports ranging from 
1921-1972. The marine weather observations incorporated in 
the heat flux calculations were of a highly diverse nature. 
The distribution of the observations was nonrandom and highly 
biased toward fair weather. Also problems existed due to 
sampling errors and methods of computation. Thus, large 
errors were possible in the heat flux estimates. Errors of 
10% in each of the components of heat exchange could 
conceivably have resulted in an error of anywhere from 
10 to 70% in the net heat flux, Q(N). 

Q(N)=Q(S) - Q(B) - Q(E) - Q(C) 
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Despite the questionable magnitude of the long term monthly 
means, however, Nelson and Husby (1981) concluded that the 
temporal and spatial consistency of the independent surface 
heat flux estimates 'indicates that the geographic patterns 
are realistic and significant’. The wind fields were of 
similar or better quality. 

The mean fields for August of Q(S), Q(B), Q(E), Q(C), 
and Q(N) (Figs. 14-18), provide insight into heat gain and 
loss within the current, and more specifically within the 
study domain. Values at the domain center (38N, 126W) , are: 
Q(S)=216.3 Watts/M** 2 
Q(B)=32.77 Watts/M**2 
Q(E)=35.92 Watts/M**2 
Q(C)=-7.24 Watts/M**2 
Q(N)=154.9 Watts/M**2 

Standard errors of the means evaluated using 
SE(Q)=SD(Q)/SQRT(N) 

where SD(Q) is the standard deviation and N is the number of 
observations (253 within the study domain region), are: 

Q(S) : Errors 3. 80 (Watts/M*“ 2 ) 

Q(B) :Error=1.41 ( Watts /M** 2 ) 

Q(E) :Error=3.28 (Watts/M**2) 

Q(C) :Error=1.47 (Watts/M’’'*2 ) 

Q(N) : Error =5. 20 (Watts/M**2) 

As expected, the incoming solar incidence corrected for 
clouds and the sea surface albedo dominates the heat 
budget. The effective back radiation and the latent heat 
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Figure 14, Long term mean August incident solar 
radiation: 1901-1972 [Nelson and 

Husby, 1983] 
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Figure 15. Long term mean August effective back 
radiation: 1901-1972 [Nelson and 

Husby, 1983] 
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Figure 16, Long term mean August sensible heat 
flux: 1901-1972 [Nelson and Husby, 

1983] 
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Figure 17. Long term mean August latent heat 

flux: 1901-1972 [Nelson and Husby, 

1983] 
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Figure 18. Long xerm mean August net heat 
exchange: 1901-1972 [Nelson 

and Husby, 19833 
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term are essentially of equal climatological value within 
the region of interest. As will be seen in the following 
discussions, the climatological value of net heat flux is 
far greater than that computed from acquired synoptic 
atmospheric data due to frequently overcast conditions. 

The climatological Q(N) is, however, somewhat less 
(approximately 30 Watts/m**2) than the value 'backed out’ 
from the transition depicted in the acquired TOPS/TOPS-EOTS 
profiles. Heat fluxes and wind stress values, the boundary 
conditions required for the model simulations, are to be 
further discussed in Chapter IV. 

Wind speed cubed, thermocline strength, and mixed layer 
depth fields are also included here to demonstrate the 
spatial and seasonal characteristics of wind generated 
turbulent kinetic energy and the associated upper ocean sta- 
bility based on Husby and Nelson (1982). Winds used in the 
computations of the wind speed cubed fields were the same as 
winds used in the calculations of wind stress and wind 
stress curl. The wind speed cubed fields were computed as an 
indicator of the turbulent kinetic energy transferred from 
the air to the upper ocean, a common procedure. 

In producing the other two fields, Husby, et.al. (1982) 
were faced with some of the typical problems one encounters 
when working with mixed layers, more specifically with 
interpreting temperature profiles (mainly XBT). It was 
necessary to establish the criterion for the extraction of 
the desired features. A typical vertical temperature profile 
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with characteristics of interest is depicted in Fig. 3. (This 
profile and some actual XBT profiles will be discussed in 
more detail in Chapter III, Section A.) That BT profiles 
comprising their data base were digitized at five meter 
intervals plays an important role in the type of criterion 
used. They chose to select their mixed layer depths as the 
upper limit of the first five meter interval in which the 
temperature gradient exceeded -0.3C/5m (-0.06C/m). The 
bottom of the thermocline was obtained by choosing the upper 
depth of the two successive five meter intervals within 
which the temperature gradient did not exceed -0.3 /5m. The 
strength of the thermocline was then obtained by taking the 
difference of the temperature at the bottom of the thermo- 
cline from the mixed layer temperature. These values were 
then summarized into one degree squares and mapped. (Fig. 19) 
The practice of using the wind speed cubed as an indica- 
tor as described above does not seem to apply to the study 
domain in August, for in using an above average estimate of 
6 m/s wind, the resulting wind speed cubed is 216 (m/s)"*3; 
a value quite a bit less than that depicted in the climatology. 
The mean observed mixed layer is, however, quite consistent 
with the corresponding climatological mixed layer depicted in 
Fig. 19. The discrepancy between turbulence indicator (wind) 
and mixed layer depth leads one to consider the probable 
importance of convective activity and horizontal advection 
within the region. 
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II. THE OBSERVATIONS 



A. DATA ACQUISITION 

The analyzed data were acquired aboard the R/V ACANIA 
in support of the Pilot Ocean Prediction Study of 
California Current Eddies, which has since become known as 
the OPTOMA program sponsored by ONR. The cruises designated 
CCSl and CCS2 (Legs I and II) covered time intervals 8 to 
12 March, 30 July to 5 August, and 9 to 14 August 1982, 
respectively. (For the convenience of this thesis, CCSl 
will be called CCSl, CCS2 legs I and II will be called CCSII 
and CCSIII, respectively.) Each had as its domain of 
interest, a region off Point Arena centered at approximately 
38N, 126W. CCSII was devoted to sampling the largest domain, 
approximately 180 km square with line spacing of ca. 44km. 
and station spacing of approximately a quarter of the 
baroclinic Rossby radius of deformation, i.e. ,9 km. This 
allowed for better resolution over a scale neither so large 
as that of the CalCOFI data nor as small as that used during 
CCSl. (120 km with line spacing of ca . 18 km and station 
spacing of 9 km. ) CCSIII was for the most part a resampling 
of the CCSII grid. Station spacing on a whole was maintained 
at small enough intervals (ca. 9 km for August cruises) to 
minimize aliasing by the smaller scale activity. 
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On all three cruises, continuous sampling was achieved 
using the shipboard DAS (Digital Data Acquisition System) . 

T-4 XBT's (expendable bathythermographs) maximum depth of 
450 m, were launched at each station. (Figs. 1, 2) 189 XBT's 

were deployed on CCSI , while approximately twice that number, 
353, were deployed in the combined latter two cruises. No 
CTD casts were made during CCSI. Interspersed along the 
track on CCSII and CCSIII were 500 and 1500m CTD casts. 

During CCSII and CCSIII, 100 CTD casts were made prior to 
equipment failure due to rough weather. In addition to 
logging these vertical profiles, DAS stored station and 
interstation "underway" values of position, time, ship's 
speed, ship's roll, and surface atmospheric data; such as, 
dewpoint temperature, air temperature, barometric pressure, 
wind direction, and wind speed. Fractional cloud cover was 
recorded in the ship's log on an hourly basis by the deck 
watchstander . 

Winds for each day of the August cruises in the vicinity 
of the cruise tracks (nine gridpoints) were supplied by FNOC. 
(Figs. 20, 21) In addition, TOPS and TOPS-EOTS produced 
ocean temperature profiles associated with eight of the nine 
points of interest which were archived by FNOC for later 
evaluation with in situ data. (The ninth point was located 
on land . ) 

Satellite IR imagery was available from NOAA/NESS 
(Redwood City) for cruise planning and preliminary analysis. 
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Figure 20. CCSII Cruise track and the nearest 
TOPS/TOPS-EOTS gridpoints 
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Figure 21. CCSIII cruise track and the nearest TOPS/ 
TOPS-EOTS gridpoints 
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Higher quality imagery was supplied by SOSF (Scripps Ocean 
Satellite Facility) for post analysis. 

B. DESCRIPTION OF DATA AVAILABLE 

1. CCS I 

The XBT data are of primary interest from the March 
cruise. 189 XBT’s were deployed, of which 180 were 
'successful*. Each trace was digitized by FNOC because 
DAS had malfunctioned. Final selection left 105 XBT’s 
within the grid region. Data files consisted of date, time, 
latitude, longitude, and vertical profiles of depth versus 
the corresponding temperature and climatologically-based 
salinity . 

2. CCSII and CCSIII 

These cruises contained more data of present interest. 
The XBT count was 257 after editing. Digitization by FNOC 
in this case was not necessary. Of these XBT ' s , 198 were 
located within the grid. 

The data, similar to that of CCSI with the exception 
that the salinity was not climatogically-based , were recorded 
on DAS cartridges. Data from these tapes were subsequently 
made accessible from CMS (Coversational Monitor System) on 
the NPS IBM 3033. The original traces and XBT logs were used 
for comparison and error checks. 

DAS log data, stored and accessed in the same way as 
the XBT data, though in the standard DAS header file format. 
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were the source of the majority of the atmospheric data 
acquired during CCSII and CCSIII. Numerous runs in various 
editing modes were made to delete blocks of zeros, erroneous 
values, and data not applicable to the grid region, (see 
Appendix A) 

Additional atmospheric data, such as cloud cover, 
were recorded on an hourly basis in the ship's log. Any 
additional data pertaining to the ocean, such as bucket 
temperatures, were logged in the CTD log or on the XBT 
traces . 

Finally, atmospheric as well as oceanic analyses in 
the vicinity of the experimental domain were archived by FNOC. 
They consisted of daily wind values, one for each of the 
nine selected FNOC gridpoints, as well as daily analyses and 
24-hour forecasts of the ocean temperature profiles for the 
eight gridpoints located over the ocean. Since TOPS/TOPS-EOTS 
runs on a Northern Hemispheric polar stereographic grid with 
a grid spacing of approximately 300 km in the midlatitudes, 
it was not possible to obtain analyses from numerous points 
in the cruise region. In fact, of the eight points located 
over the ocean, none of them fell directly within the cruise 
domain. (Figs. 20, 21) The use of the standard FNOC gridpoint 
was an unnecessary limitation that was imposed for convenience. 
(Subsequently, it was learned that FNOC has the capacity to 
interpolate their analyses to produce a vertical temperature 
profile at a specified point.) The analyses and TOPS 
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predictions used in this study were from the point closest 
to the study domain. Based on the SST climatology (Fig. 22), 
the chosen FNOC gridpoint at 39.05N, 127.75W was not 
expected to have a SST much greater (at most 1C) than that 
of the XBT’s within the cruise area. 

C. FORCING FUNCTIONS FOR THE GARV70 0D AND MELLOR LEVEL-2.5 

MODELS 

The version of the Garwood model utilized required an 
initial vertical temperature profile as well as three-hourly 
values of air temperature, dewpoint temperature, wind speed, 
and cloud cover for the computation of the boundary conditions 
to be used in both the Garwood and the Mellor 2.5 models: 
the momentum and heat fluxes. 

For the first run, initialization was done with the FNOC 
analyzed ocean temperature profile at 39. ON, 127. 75W. This 
initialization permitted comparison of the forecasted TOPS/ 
TOPS-EOTS ocean temperature profile with that obtained from 
the Garwood model over the same time span. All other runs 
were made using XBT 96 as the initial temperature profile. 

The three-hourly air temperature was obtained by 
averaging the DAS underway data. Averages were made over 6 
to 299 values. The resulting standard deviations range from 
a high of 1.67C to a low of 0.0005C. Only three values had a 
standard deviation over 1 . OC . Averaging was done under the 
assumption that the air temperature varied little over the 
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Figure 22. Average SST for last half of July 1982 
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relatively small cruise grid at any given time. In addition, 
it was assumed in the case of missing or erroneous data 
that, at a given time of day, the values did not vary much 
over the course of a couple of days. Thus, for missing 
0001-0300 readings, the 0300 input to the model was 
assumed equal to that of the previous day. (Fig. 23) 

The dewpoint temperatures were to be acquired in a 
similar manner, however, faulty instrumentation precluded 
this. Dewpoint temperatures were instead computed using a 
technique described by Bolton (1980). (Fig. 24) To use the 
method, relative humidities (also unavailable) were required. 

As recommended by Davidson (personal communication) , a 
relative humidity of 90% was assumed throughout to correspond 
with the foggy/hazy conditions which existed intermittently 
through the cruises . 

Wind speeds (Fig. 25) from the nearest FNOC point were 
used consistent with the location of initial values . The 
region of interest was considered to be too small to require 
interpolation of winds. (Fig. 26) Ship winds available from 
DAS are relative winds , but true winds could be calculated 
by correcting for ship's velocity. This was not necessary 
for the runs described. Instead, the applicability of the 
FNOC winds was checked by comparing them with the DAS winds 
recorded at CTD stations because the ship was then 'stationary', 
(see Appendix A) 



84 



CCS2 & 3 DAS AVERAGED AIR TEMP 




•H 

(T3 O 
4 h C/] 

t: c 

0 ) w o 
t:j c *H 

Ph O -P 
0-H 
a -p iH 
0) O 
^ C B 



(3 oaa) anruvnadwai hiv 



85 



Figure 23. Time series of three hourly averaged DAS 
temperatures used in computing forcing f 
the Garwood and Mellor Level-2.5 model s 
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Figure 26 . Time series of interpolated FNOC analyzed wind 
interpolated to the necessary time step 



Sea surface temperatures were obtained from individual 
XBT ' s along the track. (Fig 27) The XBT ■ s used were those 
launched within half an hour of the necessary input time. 

The final input, cloud cover in octals , was acquired 
from hourly observations, and averaged over three-hour 
intervals. (Fig. 28) 

The boundary conditions, i.e., heat and momentum fluxes, 
applied were computed within the program. (Figs. 29, 30) As 
will be discussed later, the average daily net heat flux 
forcing the Garwood and Mellor models over the thirteen days 
was +4.25 W/m**2. (where positive is upward flux, i.e., the 
ocean was losing heat to the atmosphere) 

Additional information required by the Garwood model 
includes various coefficients; such as, the extinction 
coefficient and the drag coefficient. For the most part, they 
were set from previous model usage. The Garwood model con- 
tains two 'tuning parameters' (resulting in essentially two 
degrees of freedom) which had to be addressed prior to model 
simulations with each 'new host point' [Garwood, et.al., 

1982] The parameters, AM and AZR had already been set upon 
receipt of the model. (A tuning analysis is presently 
underway within the Ocean Prediction Laboratory. The procedure 
and results will probably be documented at some future time.) 

Coefficients, such as the diffusion coefficients, used 
in the Mellor level-2.5 model are 'universal' constants, and 
so required no alterations. 
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Figure 27 . Time series Of XBT SST used to compute the 
functions for the Garwood and Mellor Level-2 
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Figure 28. Time series of cloud cover in octals 
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Figure 30. Time series of friction velocity cubed (cm/sec) 
(Julian days correspond with 31 July-14 August 
1982) 



III. STATISTICS 



A. DETERMINATION OF THE OBSERVED VARIABLE MIXED LAYER DEPTH 

Obtaining mixed layer depths from DAS or FNOC digitized 
profiles or even from the original XBT traces is no easy 
task. Criteria employed vary greatly, and there is a ten- 
dency to change one's criterion with each new exception that 
arises. Further complexity occurs with the desire to isolate 
new (shallow) mixed layer depths in addition to the numerous 
old (deep) ones. 

The new mixed layer is a layer formed over the old one 
due to the mixing being insufficient to reach the level of 
the old (previously formed) mixed layer. This may be the 
result of weak winds and surface heating. The old mixed 
layer will remain in place until turbulent diffusion of the 
interface or some other process dissipates it or strong mixing 
deepens the layer. The main reason for this is that the 
Second Law of Thermodynamics does not permit the unmixing of 
the fluid. An 'excellent' example of an old mixed layer 
recorded during CCSI is depicted in Fig. 31. (This layer is 
also the new mixed layer for it has replaced the old one 
which would have existed at an earlier time.) An old mixed 
layer extending to approximately 63m with an overlying new 
mixed layer extending to approximately 41m is depicted in 
Fig, 32. In Fig. 3, h^ would be considered a new mixed layer 
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XBT trace depicting an ’’excellent” 
new/old mixed layer 
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Figure 32. XBT trace depicting an "excellent" new 
and old mixed layer 
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overlying the old winter mixed layer, h2- New mixed layers 
do not depend solely on seasonal or synoptic changes . 

Numerous cases of diurnal mixed layers were noted wherein 
the cooling of the surface at night caused the upper water 
column to become unstable and consequently, mixed, forming a 
new mixed layer. Meanwhile, during the day, the surface is 
heated, possibly producing a warm shallow layer over the old 
mixed layer, i.e., a new thermocline is formed. 

New mixed layers may also be formed under a large number 
of other conditions. A storm, for example, will most likely 
cause a new mixed layer to form which will quite possibly 
be deeper than the old one. In this case, the new mixed 
layer supplants the old one, and becomes the old mixed layer, 
too . 

Eventhough the old mixed layers are of importance to Navy 
applications, for their role in surface acoustic ducting, the 
new mixed layers are the layers of importance in mixed layer 
modeling, i.e., that is the layer predicted and modified by 
some bulk (layered) models in use today. Thus, one must 
find a suitable method for defining both of these significant 
features, as well as sea surface temperature, depicted in an 
XBT profile. 

The first attempt to obtain mixed layer depths from the 
temperature depth profiles was made by computer comparison 
of three sequential points on the trace for a 0.2C change. 
(The value 0.2C was chosen for that is the accuracy of the 
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new XBTs; more precisely, 0.12C-0.31C [Heinmiller, 1982]. 

It is the value presently in use in producing the TOPS/ 
TOPS-EOTS contoured mixed layer depth product.) This 
technique proved to be satisfactory except in cases of a 
surface thermocline, a very weak new mixed layer, no mixed 
layer, and an old mixed layer with a noticeable gradient 
(though still possessing the characteristic ’knee’). A 
similar technique was used by Thompson (1975) with the 
additional technique of starting his sequential comparison at 
three meters depth to avoid the surface thermocline problem. 
With either procedure one would expect to obtain a depth, no 
matter whether a mixed layer existed or not. (The depth for 
the 0.2C change obtained may be in the middle of a gradient.) 
Such was not the case, however, due to the digitization 
scheme used by FNOC wherein only depths of inflection points 
were listed. Comparison of the surface temperature with the 
temperature of the first inflection point, the point at the 
base of the shallow surface thermocline, may yield a greater 
than 0.2C change. This would in turn lead to 0.0m being 
denoted the mixed layer depth when a very noticeable mixed 
layer might be apparent below the shallow thermocline or 
possibly due to the response time of the XBT recorder. 

These possible deficiencies made it necessary to return 
to the original traces and logs. Each trace was examined and 
redefined in terms of temperatures and depths of the features 
of interest, for example, surface temperature, the temperature 
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and depth of the base of the surface thermocline, the 
temperature and depth of each sequential possible mixed 
layer, and the temperature and depth of a deeper point chosen 
to define the seasonal thermocline. Final mixed layer depth 
values , as well as other values , were chosen by comparison 
of the computerized analysis, the traces, and the XBT logs. 
The computerized analysis was given the most weight, followed 
by results obtained from inspection of the traces. The XBT 
logs were given the least emphasis due to their variation 
with the experience, knowledge, and judgment of the observer 
and the not always optimal conditions under which they were 
recorded . 

B. STATISTICAL CHARACTERIZATION 

The several variables derived from the XBT profiles by 
the techniques previously mentioned were used as a basis for 
a statistical characterization of each cruise. Basic pro- 
perties such as the mean, minimum, maximum, median, and 
standard deviation were examined. They are summarized in 
Table III. Two types of correlations coefficients were 
calculated, as well as Student-t confidence intervals, 
regression equations, and histograms and maps of SST, bucket 
temperature, old mixed layer depth (OMLD) , old mixed layer 
depth temperature (OMLDT), new mixed layer depth (NMLD) and 
new mixed layer depth temperature (NMLDT). Only a limited 
selection of statistical quantities will be presented to 
convey the character of the upper ocean during the various 
cruises . 
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The first variables to be discussed are the SST's for 
the three cruises and the bucket temperatures from the first 
cruise. A total of 105 XBT ' s were available from the CCSI 
grid for the statistical characterization. Both XBT SST and 
bucket temperature were recorded with a comparison revealing 
the XBT SST field was overall slightly warmer than the bucket 
temperature field, i.e., O.^C. (Fig. 33, Table III) A 
regression equation or least squares linear equation of the 
bucket temperature with respect to the XBT SST was calculated 
to be: Y=.678 + .926(X1) where X1=XBT SST. The R squared 

value, a measure of how well the regression equation fits 
the data, was 79.6%. (100% is a perfect fit. The R squared 
value is defined by (100) (sum of squares due to regression/ 
total sum of squares . ) ) 

The two correlation coefficients computed were the Pearson 
product moment coefficient and the autocorrelation function 
with respect to distance. The former calculated for XBT SST's 
versus bucket temperatures, resulted in a large correlation 
coefficient, .892., though not 'large enough' considering 
the fact that the values should be exactly the same (except 
for measurement noise). The spatial autocorrelation function 
was computed using binned XBT pairs. Out of 210 km, the 
radial resolution for XBT pairs within a bin was 15 km. 

Beyond 210 km, it was increased to 25 km. XBT SST for CCSI 
(Fig. 34) had a correlation distance of approximately 30 km, 
using a 0.2 correlation as a reference. The correlation 
distance for the bucket temperature was slightly less, i.e., 

27 km. (Fig. 35) 
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Figure 34. Autocorrelations function for CCSI bucket temperature 
with respect to distance; results depicted are for 
detrended data; also indicated are the number of data 
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Figure 35. Autocorrelation function for CCSI XBT SST with res 
distance: results depicted are for detrended data 

'raw' data (with trend); also indicated are the nu 
data pairs 



TABLE III 



BASIC STATISTICS 





CCSI 


CCSII 


CCSIII 




SST Bkt 


SST 


SST 


Mean 


11.5 11.3 


15.9 


16.3 


Stddev 


.302 .314 


. 91 


.703 


Min 


10.9 10.4 


12.9 


13.9 


Max 


12.4 12.0 


18.0 


17.1 


Median 


11.5 11.3 


15.8 


16.5 


T 90% 


(11.4,11.5) (11.3,11.4) 


(15.7,16.1) 


(16.1,16.4) 


C.I.* 


NMLDT 


NMLDT 


NMLDT 


Mean 


11.3 


15.4 


15.9 


Stddev 


.324 


.822 


.744 


Min 


10.8 


12.3 


11.9 


Max 


12.2 


16.7 


16.9 


Median 


11.3 


15.6 


16.2 


T 90% 


(11.3,11.4) 


(15.3,15.6) 


(15.8,16.1) 


C . I . “ 


NMLD 


NMLD 


NMLD 


Mean 


39.1 


15.1 


16.3 


Stddev 


25.5 


8.3 


8.1 


Min 


1.0 


4.0 


5.0 


Max 


82.0 


55.0 


50.0 


Median 


33.0 


13.0 


15.0 


T 90% 


(34.9,43.4) 


(13.6,16.6) 


(14.9,17.7) 


C.I. * 


OMLDT 


OMLDT 


OMLDT 


Mean 


11.2 


14.7 


15.5 


Stddev 


.287 


.927 


.801 


Min 


10.5 


11.8 


11.9 


Max 


11.7 


16.4 


16.6 


Median 


11.2 


14.8 


15.7 


T 90% 


(11.1,11.2) 


(14.6,14.9) 


(15.4,15.6) 


C.I.* 


OMLD 


OMLD 


OMLD 


Mean 


64.8 


37.8 


29.5 


Stddev 


12.0 


13.4 


12.2 


Min 


27.0 


6.0 


8.0 


Max 


95.0 


69.0 


65.0 


Median 


63.5 


38.0 


28.0 


T 90% 


(62.8,66.7) 


( 35 . 5,40.1) 


(27.5,31.6) 


C.I.* 

Unites 


in table: all temperatures: degrees Celsius 


* T 9 0% 


all depths : meters 

C.I.: T Interval With 90% Confidence 
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During the latter two cruises, bucket temperatures v/ere 
not recorded routinely. The CCSII data set consisted of 97 
XBT's. As was to be expected, the SST values were much 
higher than those recorded for the March cruise. In fact, 
the minimum CCSII temperature, 12. 9C, exceeded the maximum 
CCSI , 12.4C. The range of temperatures for CCSII, and thus 
the standard deviation in this case was quite large. A 
reason for the significant range, 12.9 - 18. OC, was the 
presence of cold eddies and cool anomalies near the peak 
of summertime warming. 

The autocorrelation functions for SST for CCSII were 
plotted with and without a linear trend removed. Fig. 36. 

This was done to get some idea of the existence and effects 
of a trend in the data. Also shown is a 95% confidence 
region computed by +/- sqrt(2)/N where N is the number of 
paris in the bin. (Using this N instead of the total N 
results in a more conservative confidence region.) Results 
indicate a slight increase of the correlation distance 
compared to the March cruise, i.e., 32 km for the detrended 
data and 38 km without the trend removed. 

The mean SST for CSSIII was still higher, 16. 2C, as might 
be expected for later in August. The range, however, fit 
within that of the second cruise. This may be related to 
the fact that CCSIII did not cover as large a region as CCSII. 
The correlation distance for this cruise was 40 km for the 
detrended data and 43 km without the trend removed. (Fig. 37) 
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This might indicate a more homogeneous area as compared to 
the larger grid transited during CCSII. 

The new mixed layers for the March cruise ranged from 
0.0 m (when no new mixed layer was observed) to 82 m. (Table 
III) The largest Pearson products moment correlation obtained 
was for the NMLD’s and OMLD ' s • This resulted due to a number 
of cases when the NMLDw;as identical to the OMLD, i.e., in 
cases of intense turbulence. 

For CCSII, 84 XBT ’ s possessed a NMLD . The NMLD’s were, 
however, quite small, approximately half of those noted in 
CCSI , (Fig. 38). CCSIII had 94 instances of NMLD’s. The 
distribution in the histogram appears to be Raleigh in 
nature . 

Histograms of the OMLD’s for the cruises are shown in 
Fig. 39. Most noticeable is the shift toward shallower 
mixed layers with time. In fact, the mean depth for the March 
cruise is a factor of two greater than the CCSIII mean OMLD. 
The Pearson product moment correlation coefficients for the 
OMLD in association with the majority of the other variables 
showed no dominant results. The largest coefficient obtained 
was with respect to the NMLD*, the probable reason for this 
was already discussed. For CCSI, the correlation distance 
was merely 15 km. (Fig. 40) For CCSII, both the data with 
and without the trend removed, have a correlation distance 
of only approximately 14 km. (Fig. 41) CCSIII possessed old 
mixed layers that were more highly correlated than those 
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Figure 39. 



Histograms Of OMLD (meters) for CCSI, 
CCSII , and CCSIII . 
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Figure 40. Autocorrelation function for CCSI old mixed 
layer depth 
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Figure 41. Autocorrelation function for CCSII old mixed 
layer depth 



previously examined. (Fig. 42) The correlation distance for 
the detrended data was approximately 35 km; without the 
trend removed the correlation distance was approximately 50 
km. The increase in the correlation distance from CCSII to 
CCSIII may be due to the fact that the former cruise covered 
a wider and more diversified domain than did the latter. 

Contours of OMLD are included for CCSII and CCSIII to 
give an idea of the variability of mixed layers (a quasi- 
synoptic depiction) , and to bring to light a common problem 
encountered in trying to contour ’real’ data. The data must 
be placed on a regularly spaced grid prior to using contouring 
graphics packages. One would think the answer to the problem 
would be to merely interpolate the data to the new grid. 
Problems arise, however, from data set to data set depending 
on the noise or gradients in the data. Certain schemes seem 
to work well for some sets and not to work for others. The 
enclosed computer plotted contoured fields are the result of 
interpolation by trying to fit a five degree polynomial to the 
data. (Fig. 43) The results do not have much worth as seen 
when compared with the hand contours. (Fig. 44) Only the 
locations of some of the minimum/maximum centers are retained. 
In the future, the technique of optimal interpolation should 
be applied to the data for improved interpolated values and 
consequently, more appropriately smoothed contoured fields. 
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Figure 42. Autocorrelation function for CCSTII old mixed 
layer depth 
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Figure 43. Computer ciontoured old mixed layer 
depths from CCSIII data were 
- interpolated to a regular rotated 
grid and contoured; rotation origin: 
37.33N, 126. 25W; subplot origin: 

38. ION, 126. 38W 
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Figure 44. Hand contoured old mixed layer depths 
(meters) from CCSIII 
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IV. ANALYSIS OF OBSERVED AND SIMULATED 
TEMPERATURE STRUCTURE 



The first run was of the Garwood model initialized with 
the TOPS predicted profile for 1 August. The initial mixed 
layer depth was set at 20m, and the initial SST chosen was 
17. 7C (the SST of the TOPS analysis profile). The 24 daily 
profiles were averaged so as to produce one 'representative' 
profile per day. The average profiles along with their 
standard deviations compose Figs. 45-58. The standard 
deviation was nonzero only in the region at and/or above the 
base of the mixed layer, because the Garwood model does not 
alter regions below the furthest extent of the mixed layer 
base. For the 13-day cruise period, the Garwood model 
cooled the mixed layer by approximately two thirds of a 
degree, and deepened the mixed layer to approximately 37m. 

A substantial amount of entrainment occurred from 2 to 3 
August, with the simultaneous increase in wind speed (7.8 m/s), 
and decrease in the SST at OOZ. Thus, there was increased 
forced convective activity due to increased wind stirring and 
increased free convection due to surface ccoling. Significant 
mixing also occurred over the time spans; 7 to 8 , 9 to 10, 

12 to 13 August with some entrainment also occurring from 
11 to 12 and 10 to 11 August, but not as pronounced as that 
during the previously stated times. Not much occurred during 
the other time spans. 
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Figure 45. Overall average profile from Garwood 

model (Temperature in degrees Celsius) 
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AVE MODEL PROFILE, AUG 1 



TEMPERATURE 
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Figure 46. 



Average 1 August profile from Garwood 
model (Temperature in degrees Celsius) 
Initialized with FNOC analyzed profile 
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AVE MODEL PROFILE, AUG 2 



TEMPERRTURE 
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Figure 47. Average 2 August profile from Garwood 
model, Initialized with FNOC analyzed 
profile (Temperature in degrees Celsius) 
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AVE MODEL PROFILE, AUG 3 



TEMPERRTURE 
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Figure 48. Average 3 August profile from Garwood 
model Initialized with FNOC analyzed 
profile (Temperature in degrees Celsius) 
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AVE MODEL PROFILE, AUG 4 



TEMPERflTURE 
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Figure 49. Average 4 August profile from Garwood 
model. initialized with FNOC analyzed 
profile (Temperature in degrees Celsius) 
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AVE MODEL PROFILE, AUG 5 



TEMPERRTURE 



0.0 5.0 10.0 15.0 20.0 25.0 




Figure 50. Average 5 August profile from Garwood 
model. Initialized with FNOC analyzed 
profile (Temperature in degrees Celsius) 
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AVE MODEL PROFILE, AUG 6 



TEMPERATURE 
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Figure 51. Average 6 August profile from Garwood 
model. Initialized with FNOC analyzed 
profile (Temperature in degrees Celsius) 



124 



AVE MODEL PROFILE, AUG 7 



TEMPERRTURE 



0.0 5.0 10.0 15.0 20.0 25.0 




Figure 52. Average 7 August profile from Garwood 
model, Initialized with FNOC analyzed 
profile (Temperature in degrees Celsius) 
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AVE MODEL PROFILE, AUG 8 



TEMPERATURE 



0.0 5.0 10.0 15.0 20.0 25.0 




Figure 53, Average 8 August profile from Garwood 
model. Initialized with FNOC analyzed 
profile (Temperature in degrees Celsius) 
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AVE MODEL PROFILE, AUG 9 



TEMPERATURE 
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Figure 54. Average 9 August profile from Garwood 
model Initialized with FNOC analyzed 
profile (Temperature in degrees Celsius) 
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AVE MODEL PROFILE, AUG 10 



TEMPERATURE 



0.0 S.O 10.0 IS.O 20.0 2S.0 




Figure 55. Average 10 August profile from Garwood 
model. initialized with FNOC analyzed 
profile CFemperature in degrees Celsius) 
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Figure 56,. Average 11 August profile from Garwood 
model , Initialized with FNOC analyzed 
profile (Temperature in degrees Celsius) 
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AVE MODEL PROFILE, AUG 12 



TEMPERRTURE 
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Figure 57. Average 12 August profile from Garwood 
model, initialized with FNOC analyzed 
profile (Temperature in degrees Celsius) 
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AVE MODEL PROFILE, AUG 13 



TEMPERfiTURE 
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Figure 58. Average 13 August profile from Garwood 
model. Initialized with FNOC analyzed 
profile (Temperature in degrees Celsius) 
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The time-depth contours of temperature in Fig, 59 are 
another way of indicating the deepening and shoaling of the 
mixed layer depth. It is quite obvious that the depths 
below the mixed layer are not affected by the model. 

The TOPS profiles, on the otherhand, are affected at all 
levels (down to 500m). The profiles are more complex and 
harder to interpret, because of the variations occurring 
over the entire profile. The process of energy addition/ 
subtraction to the mixed layer is no longer as obvious as 
with a bulk model. Overall, changes over the 13 days, 
however, clearly show a substantial increase in the mixed 
layer temperature; the SST increased approximately one 
degree, while the mixed layer depth increased to approximately 
40m. (Figs. 60, 61) (TOPS and TEOTS profiles for the other 
days may be found on the composite figures to be discussed 
in the last section.) 

A comparison of the contours for the TOPS predicted 
profiles (Fig. 62) with those of the Garwood model gives 
further indications of the difference in results of the 
diffusion and the 'slab' or integral models. For these 
contours to resemble those obtained from the Garwood model, 
the mixed layers would indeed have to be 'slab-like'. (The 
'plateaus' located in the center of the contour field are due 
to repeated profiles, i.e., profiles for 6 and 7 August were 
not obtained from FNOC.) 
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Figure 59. Contoured Garwood model simulated T/Z Field 



FNOC TOPS&TEOTS PROFILES 

39. 1N,1 27.7 W:D ATE-82080 1.00Z 

TEMPERATURE 




Figure 60. TOPS and TEOTS profiles for 1 August 

1982. (Temperature in degrees Celsius) 
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FNOC TOPS&TEOTS PROFILES 
39. 1 N, 1 27.7W:DATE-8208 1 3.00Z 

TEMPERATURE 



0.0 S.O 10.0 15.0 20.0 25.0 




Figure 61. TOPS and TEOTS profiles for 13 August 

1982 (Temperature in degrees Celsius) 
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Figure 62. Contoured T/Z field composed of TOPS profiles 



Both of these sets of profiles and contours were compared 
with those obtained from daily averaged XBT profiles . The 
variability of the upper water column as noted in the latter 
two sets of profiles was not as apparent in the series of 
averaged XBT profiles, as expected. The trend of the latter 
XBT ' s, averaged in space and time, had a tendency to wander 

quite a bit. For example, the SST from 1 August to 2 August 

jumped 1C with a corresponding increase of 0 . 3C at 400m. 
(Figs. 63, 64) Similar changes occurred for other profiles. 
For example, the 13 August profile indicated a -1.3C change 
in SST and a -0.7C change in the temperature at 400m since 
12 August. (Figs. 65, 66) (Average XBT profiles for the 

other days may be found on the composite figures to be 

discussed in the last section.) Mixed layers and their tem- 
poral variations are not well represented by these profiles 
due to the spatial heterogeneity within the study domain. 

It takes a stretch of the imagination to identify a mixed 
layer in som.e of the profiles. Once identified, the mixed 
layers do not tend to exhibit the characteristics expected 
for mixed layer evolution, e.g., the deeper layers are not 
necessarily cooler. The characteristics seem to be more of 
a 'random' nature, i.e. they are affected by factors other 
than those of mixing and heating. 

The use of average XBT profiles is somewhat of a contro- 
versial issue. One side holds that the averaging of XBT ' s 
has some worth in that it will cancel some of the spatial 
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AVE XBT PROFILE, AUG 1 



TEMPERATURE 
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Figure 63. Average XBT and standard deviation plot 
for 1 August 1982. (Temperature in 
degrees Celsius) 
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AVE XBT PROFILE, AUG 2 



TEMPERATURE 

0.0 5.0 10.0 15.0 20.0 25.0 




Figure 64. Average XBT and standard deviation plot 
for 2 August 1982. (Temperature in 
degrees Celsius) 
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AVE XBT PROFILE, AUG 12 
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Figure 65. 



.Average XBT and standard deviation plot 
for 12 August 198 2 . (Temperature in 
degrees Celsius) 
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and temporal effects. The other contends that the averaging 
of XBT ’ s destroys the vertical structure, and, therefore, 
should not be the method by which one obtains a single daily 
profile with which to compare the model simulations. The 
latter believes that a median XBT (based on various cri- 
teria) would be a much better ’measuring stick’ for model 
simulation success . Having gone through the abundant supply 
of XBT traces, comparing the averaged profiles with the 
series, the author is convinced of the latter viewpoint, 
though not contradicting that of the former. The corollary, 
then, is that the thermal structure is not robust over 
modest distances of a day’s transit in a slow ship (ca. Sms ^). 
It is also suggested that greater care should be taken in 
grouping the XBT profiles in space-time for analysis. 

Indicated on the righthand panels in Figs. 67 and 68 are a 
series 'of hypothetical XBT traces for two cases: one 

containing XBT ’ s with pronounced jumps at the mixed layer 
bases, and the other containing XBT ’ s with gradually sloping 
mixed layers. The lefthand panel of each figure depicts the 
average XBT trace for the series, as well as a median XBT. 

(The average profile is the shallower of the two.) From 
these it can be seen that averaging distorts the vertical 
structure, making the mixed layer depth shallower, and the 
temperature at the base of the mixed layer colder. 

Median XBT profiles were almost as hard to define as 
were the mixed layer depths, especially over a substantial 
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Figure 67. Averaging XBT's with steep thermoclines A) Right panel: 
XBT's to be averaged; B) Left panel: Average 

(indicated by arrow) and median XPTts 
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Figure 68. Averaging XBT's with less steep thermoclines A) Right 
panel: XBT’s to be averaged; B) Left panel: Average 

(indicated by arrow) and median XBT’s 



time and distance. XBT traces for a day varied from having 
’excellent bulklike’ mixed layers to having no mixed layers. 
The SST’s did not vary in any predictable pattern, i.e., 

SST’s of XBT traces without mixed layers were in fact colder 
than the SST’s of the traces with mixed layers. The transit 
of the ship across numerous fronts, jets, and eddies would 
be one explanation for these anomalies and ambiguities. 

Instead of obtaining a daily median profile, profiles 
for the same time each day were used as the ’observed’ 
reality. The time chosen was 1800Z or 1000 local time to 
allow for the largest sample of XBT ’ s , and to avoid periods 
of major diurnal heating or cooling. Fig. 69 is the T/Z 
contour for the XBT ’ s closest to 1800Z each day, approximately 
within 30 minutes. There were no XBT’s available for the 
designated time on 5, 6, 7, or 8 August. Therefore, the 
profiles for 5 , 6 August were made identical to that of 4 
August. The profiles for 7, 8 August were made identical to 
the profile for 9 August. The ’wavelike’ nature and the two 
events present in the average XBT T/Z contour field at 2, 12 
August are also located in the 1800Z XBT T/Z contour field. 
Associated with these peaks or troughs, one might expect to 
find a change in the mixed layer due to fluctuations in the 
whole water column. Such in fact was the case for mixed 
layers of 2 and 12 August. 
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Figure 69. Contoured cemperature/depth Of 1800Z XBT 



Returning to the results of the Garwood and TOPS model 
runs discussed thus far, one must keep in mind a number of 
important points. Recall that the Garwood model simulation 
shown was initialized with a profile from outside the study- 
region and with inputs of varying degrees of accuracy. It 
was not known exactly how good the boundary/ conditions were, 
nor how sensitive that model would be to them. A sensitivity 
test was performed to decide whether the Garwood model was 
deepening the layer due to the 19.5m FNOC winds being too 
strong. This was done by decreasing the drag coefficient 
by 10%, equivalent to decreasing the FNOC winds to 95% of 
their original magnitude. As can be seen from profiles for 
1 and 13 August 1982 (Figs. 70, 71), the effects are not very 
substantial . 

In addition, it should be remembered that the TOPS model 
was initialized and run at a point outside of the cruise 
grid. It is not known what values TOPS received as initial 
or boundary conditions, only that it received its atmospheric 
forcing from the atmospheric primitive equation model and 
NOGAPS, and that eight XBT’s were used in the TEOTS analysis. 
Finally a slightly warmer SST is expected for the FNOC 
gridpoint, because it is further offshore than the study 
domain. According to the previous fortnightly mean, 
however, the difference should only be about 1C. 

To make simulations more consistent within the domain, a 
second run was performed using XBT 96 to initialize the 
Garwood model. This XBT, deployed on 4 August, was chosen 
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AVE PROFILE 
39.1N 127.7 W DATE= 8/1 

temperature: 




Figure 70. Daily average and average tuned Garwood 
profiles: 1 August (Temperature in 

degrees Celsius) 
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AVE PROFILE 
39.1N 127.7W DATE= 8/l3 



_ . TEMPERRTURE 
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Figure 71. Daily average and average tuned Garwood 
profiles: 13 August. (Temperature in 

degrees Celsius) 
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because its location was resampled on 10 August with the 
launch of XBT 216. (Fig. 72) (It was hoped that by choosing 
XBT 96 to initialize the Garwood model, its mixed layer 
evolution could be compared and evaluated with the data 
obtained from XBT 216. The initial SST used was that of 
XBT 96. An initial mixed layer depth of 20m v;as "assumed".) 

In addition, to avoid comparing simulations of two 
differently forced models, the Mellor level-2.5 model was 
obtained and modified to produce simulations using forcing 
identical to that used in the Garwood model. Profiles 
resulting from the model runs (Figs. 73-85) along with the 
applicable 1800Z XBT profiles, TOPS predicted profiles, and 
daily averaged profiles are examined. (Motes: A) repeated 

TOPS, 1800Z, and average XBT profiles during in-port period; 
B) TOPS profiles correspond to 1600 local time, i.e., 
afternoon heating exists; C) Mellor profiles are not daily 
averages . ) 

Based on the overall change in the profiles , Garwood and 
Mellor profiles cooled and deepened the mixed layer. The 
profiles seemed to be quite consistent with the observed 
data by 4 August and even more so by 10 August. The cooling 
trend in the XBT ' s seemed to be consistent with that apparent 
from XBT 96 to XBT 216. However, caution is required, 
because the role of the horizontal and vertical advection 
has not been assessed. 
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Figure 72. XBT ’ s 96 and 216 (Temperature in 

degrees Celsius) launched at approximately 
the same location on 4 and 10 August, 
respectively 
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Figure 73. August 1 composite of simulated and 
observed profiles 
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DAILY PROFILE, 8/2 
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Figure 74. August 2 composite of Simulated and 
observed profiles 
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DAILY PROFILE, 8/4 
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Figure 76. August 4 composite of simulated and 
observed profiles 
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DAILY PROFILE, 8/5 
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Figure 77. August 5 composite of simulated and 
observed profiles 
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DAILY PROFILE, 8/6 
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Figure 78. August 6 composite of simulated and 
observed profiles.. 
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DAILY PROFILE, 8/7 
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Figure 79. August 7 composite of simulated and 
observed profiles. 
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Figure 80. August 8 composite of simulated and 
observed profiles 
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Figure 81. August 9 composite of simulated and 
observed profiles 
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Figure 82. August 10 composite of simulated and 
observed profiles 
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Figure 83. August 11 composite of simulated and 
Observed profiles . 
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DAILY PROFILE, 8/12 



TEMPERATURE (C) 




Figure 84, August 12 composite Of simulated and 
observed profiles. 
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DAILY PROFILE, 8/13 
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Figure 85. August 13 composite Of simulated and 
observed profiles,. 
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The TOPS profiles, though not expected to resolve the 
small scale veriability of the boundary current region, on 
the otherhand, started out warmer than the other profiles and 
indicated additional v/arming with time. The fact that the 
gridpoint used was further offshore than the cruise region 
caused one to expect a slightly greater surface temperature 
for the FNOC profiles. That the forcing averaged 184 Watts/ 
m"*2 per day was, however, not anticipated. This heating 
is, however, believed to be part of a larger problem, viz., 
a heat bias in the atmospheric PE model. Heats biases of 
varying signs have been previously identified by Elsberry, 
et.al. (1979), Budd (1980), and Warn-Varnas , et.al. (1982) 
and shown to possess temporal and spatial scales of a month 
and tens of thousands of kilometers, respectively. Warn- 
Varnas, et.al. (1982), for example, found that a limited- 
area version of TOPS initialized and forced with FNOC fields 
revealed considerable skill during 60-day runs in comparison 
to persistence if a bias of about 150 ly/day was removed 
from the daily net surface heat flux of the domain. This 
correction proved to yield consistent results with the heat 
fluxes computed directly using atmospheric observations 
acquired by ships-of-opportunity . [Clancy, et.al., 1980] 

In another study, comparison of the FNOC heat fluxes with 
those calculated by Nate Clark indicated a cold bias believed 
to be mainly the result of excessive latent heat loss. The 
problem of heat bias is extremely evident in such data sparse 



areas as the high latitudes. For example, TEOTS SST's in 
the Sea of Okhotsk, the Bering Sea, and the Labrador Sea 
have in the past exhibited a warm bias, running as much as 
6C above those of conventional EOTS (essentially 
representative of climatology within these areas). [Clancy, 
et.al, 1982] 

In the simulations presented, the initial conditions and 
forcing for the first two days were supplied to TOPS by the 
PE model. On 3 August, NOGAPS replaced the PE model. 

However, any spurious heat available in the PE model entered 
into the simulation prior to bringing NOGAPS online. Once 
introduced into the run, it is unlikely that the heat could 
be removed without much mixing or surface cooling (due to 
the Second Law of Thermodynamics). NOGAPS’ contribution to 
the spurious heating is unknown, though it's clear that the 
heating trend was not corrected over the short run time. 

(It is as yet unknown whether NOGAPS has a heat bias of any 
kind . ) 

Examining the TEOTS SST anomalies (TEOTS - climatology) 
for some of the cruise days as well as for a number of days 
following the cruise, and following the 'return to climatology' 
adjustment installed in TOPS on 20 September (Table IV) 
revealed an excess of 2-3C in the study domain as compared 
to conventional EOTS. 
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TABLE IV 



TEMPERATURE ANOMALIES 



DATE 


TEMPERATURE ANOMALY CO 


820730 


TEOTS EOTS 

2 


820803 


NOGAPS Online 


820804 


2 -1 


820808 


2 0 


820811 


3 


820814 


2 


820819 


3 


820827 


2 


820828 


System Restart 


820910 


1.7 


820920* 


1.4 


820923 


1.9 



^Return to climatology adjustment installed 

The warmer SST's along with the an identified tendency 
for the turbulent parameter scheme used in TOPS to underpredict 
the MLD's [Martin, 1982] might be important in explaining the 
trend in the TOPS’ profiles. The difference between various 
model profiles with depth is believed to be the result of the 
strong connection TEOTS has with climatology. 
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To make a strict comparison of model results (assuming 
ideally that the model forcing is the same, as it should be), 
one should interpolate between FNOC gridpoints to the 
locations of XBT 96. 

It is felt , however, that the heating difference is not 
totally a function of distance. Erroneous heat fluxes, 
strong dependence upon climatology, and advection (associated 
with mesoscale variability) are thought to play important 
roles in the results. 
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V. SUMMARY AND CONCLUSIONS 



The Garwood and Mellor level-2.5 simulations resulted in 
cooling and deepening of the mixed layer over time, and 
exhibited a likeness to the observations. The TOPS model, 
receiving different forcing, produced warmer profiles with 
a substantial warming trend over the 13 days. 

The exhibited differences of the TOPS profiles with 
respect to the profiles produced by the other models are 
believed to be more than a result of location. They are 
believed to be due to a) a heat bias in FNOC's atmospheric 
PE model, b) too much dependence upon climatology, and 
c) advection and the subgrid scale variability existing in 
the region. 

The autocorrelation functions for SST and MLD indicated 
j a maximum correlation distance of 35 km. However, TOPS 
uses a standard FNOC grid with spacing 381 km at 60N. 

Hence, though interpolation should be carried out when looking 
at regions between the gridpoints, it seems highly unlikely 
to be able to recover the small scale variability dominating 
the region. 
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VI. RECOMMENDATIONS 



First of all, numerous recommendations are in order 
concerning the data acquisition on future cruises. Most 
importantly, all easily observed data should be recorded 
regardless of the immediate mission of the cruise, (e.g., 
wet bulb temperature) The equipment should be in a working 
status, and be sufficiently backed-up with manual logs. 
Underway analysis should be carried out in a consistent 
manner wherever possible to avoid excessive reworking of the 
data ashore. A thorough scientific crew indoctrination 
period is, therefore, in order. 

From the modeling aspect, many m.ore evaluations of TOPS 
are necessary under numerous varied conditions prior to its 
full acceptance by the scientific ' community . It is important 
that the TOPS products be interpolated to the observed 
position and tLme of the study domain prior to any further 
evaluations of the model's performance. Design tests should 
be carried out with control over the variables. Independent, 
direct measures of atmospheric forcing should be made to 
drive the research and operational models. CSimulations 
should be produced using both sets of forcing values). 
Additionally, it would be of interest to attempt to predict 
SST and MED from one realization to the next, and to examine 
the changes in the profiles in the (z,x,t) domain instead of 



170 



merely the (z,day) domain. All of these tests and comparisons 
should also be performed on the 60 km and 20 km grids of the 
future . 

In addition, other operational models and their products 
are becoming more important both in the Navy and in the 
civilian community. With their growing availability, so too 
must grow the amount and quality of verification testing 
accorded them, for they are an important asset of the future. 
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APPENDIX A 



A. DATA PROCESSING PROCEDURES 

Working with a new data set is rarely an easy task. The 
data must be edited and processed into a usable format. 

Here, the editing procedures used to produce the five 
required inputs for the version of the Garwood model at hand 
are described in more detail to better assess the model 
inputs and results . 

The procedure applied to the DAS data was as follows. 
First, the data on magnetic tape was dumped onto mass storage 
and filtered for blocks of zeros and gross errors in lati- 
tude and longitude. Following this, the data from the 
entrance and exit tracks to the cruise grid were deleted so 
as to confine the data to the study domain. Finally, the 
finer errors in temperature, wind speed, etc. were deleted. 

The resulting air temperature data were then averaged over 
three-hour intervals to yield the required eight daily air 
temperatures. Averages were computed over bins containing 
from 6 to 299 values with the maximum standard deviation for 
any resultant value being 1.67C and the minimum being 
0.0005C. Only three of the averaged values possessed 
standard deviations over 1C. Procedures used to obtain 
missing values within cruise days were based on the 
assumptions described previously. For days in port, the 
air temperatures were held at a representative constant value. 
(This could be improved upon with reasonable means to predict 
gaps . ) 
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The dewpoint temperatures were erroneous due to equipment 
failure . The required eight values of dewpoint per day were 
calculated by assuming a relative humidity of 90% and using 
the procedure described in Bolton (1980). This value for 
relative humidity was chosen to correspond with the 
intermittently hazy/foggy conditions. 

Cloud cover indications from the ship's log were either 
in tenth of sky cover or in descriptive terms. The first 
step taken was to convert the descriptors to tenths of sky 
cover using Bowditch’s American Practical Navigator and the 
experience of the R/V ACANIA's first mate. The tenths of 
sky cover were then converted to octals, and averaged 
over three hour intervals. Missing values during cruise 
days were not a problem as they had been with the air and 
dewpoint temperatures . The cloud cover during the in-port 
days was held at a constant five tenths. 

Winds could have been obtained in two ways . The DAS 
recorded winds could have been converted to true winds and 
averaged as was done with the air temperatures . The other 
way winds could have been obtained was by interpolating the 
daily analyzed TEOTS 39.05N, 127. 75W wind to the necessary 
three-hour interval. The TEOTS analyzed winds over both 
CCSII and CCSIII are depicted in Figs. 25 and 26 with the 
latter time series containing the interpolated eight three- 
hourly values. (The symbol denotes the value at 00 Zulu) 

The latter method was chosen for convenience as well as to 
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facilitate the comparison of the results of the TOPS, Mellor 
level-2.5, and Garwood models. Some of the 'CTD station 
winds’ (already assumed to be true winds due to the ’zero' 
speed of advance of the ship) were plotted against the FNOC 
archived winds to observe whether or not the latter could 
’reasonably' be applied to the cruise area. The CTD winds 
were found mainly to fluctuate around the FMOC winds except 
for the first day and a half when the 'CTD winds’ were 
exceeded by the FNOC winds. The FNOC winds were (Fig. 86) 
employed in subsequent runs. (The triangles represent 'CTD 
winds’ recorded along with suspect dataj leading to a 
possible suspicion in the winds.) A sensitivity analysis was 
done to assess the importance of the wind stress as a bound- 
ary condition, and in so doing, ascertain the magnitude of 
possible problems introduced by using winds recorded at 
19.5m vice 10.0m. The procedure and results were described 
in Chapter IV. 

The final input required was SST. The SST’s were avail- 
able from a number of sources. DAS recorded boom temperatures 
could have been averaged. This route was not chosen, because 
of the possible problems posed by the boom varying in its 
depth in the water (interaction with the ship’s wake). 
Alternatives included the use of the bucket or XBT surface 
temperatures, averaged or singly, at the required times. 

For the initial runs, the XBT SST’s recorded by DAS 
were used: this was believed to be the least subjective 
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Figure 86, Time series of FNOC analyzed winds Vs DAS logged winds 



approach. The bucket temperatures were subject to the 
experience of the watchstander , as were the XBT SST's 
recorded in the XBT logs. The validity of the latter were 
found to strongly depend upon the existence of a surface 
thermocline, either ’real’ or possibly due to the initially 
nonlinear slow response rate of the XBT recorder. The XBT 
temperatures were taken from those launched closest to the 
required input time. (usually within 30 minutes--sometimes 
linear interpolations were performed.) Averages were not 
made for the three-hour transit of the ship across possible 
fronts . 

Finally it was necessary to input an initial temperature 
profile into Garwood and Mellor level-2.5 models. Temperatures 
were interpolated to intervals of one meter to a depth of 
200M. The first runs were made using an interpolated TEOTS 
profile corresponding to 1 August. (The original TEOTS 
analyzed temperature profiles consisted of 13 temperature 
depth values. Surface values were not available; the tempera- 
ture at one meter was used at the surface.) Subsequent runs 
were made using XBT 96 as the initial temperature profile for 
the Garwood and Mellor level-2.5 models. This XBT was chosen, 
because of its proximity to the FNOC gridpoint of interest. 

More importantly, however, it was chosen because of its common 
position with an XBT, 216, launched six days later. XBT ’ s 96 
and 216 were launched on 4 and 10 August, respectively, thus 
spanning a part of both cruises. 
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APPENDIX B 



TABLE V 

VERTICAL RESOLUTION OF TOPS 
VERTICAL GRID SPACING USED IN TOPS 



Level 


DepthCm) 


1 


2.5 


2 


7o5 


3 


12.5 


4 


17.5 


5 


25.0 


6 


32.5 


7 


40.0 


8 


50.0 


9 


62.5 


10 


75.0 


11 


100.0 


12 


125.0 


13 


150.0 


14 


200.0 


15 


300.0 


16 


400.0 


17 


500 . 0 



TOPS is run on the Cyber 205 computer on a standard 
63x63 FNOC global grid. Values of T, S, U, and V are defined 
with this horizontal resolution at each of the seventeen 
levels listed above. All spatial derivatives are centered 
in space. The horizontal advection components are defined 
on a staggered grid; they are forward-differenced in time. 
Vertical eddy fluxes and W are defined midway between these 
levels. The former are differenced backv/ard in time, while 
the latter (as well as the Coriolis term) is differenced 
trapezoidally in time. 
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r 








Additionally, TOPS is quasi-three-dimensional model, 
(though a nonadvective version is kept on the Cyber 175 as a 
backup.) Advection, as previously mentioned, is by 
climatologically averaged geostrophic currents and 
instantaneous wind drift. The geostrophic currents are 
updated monthly using density fields computed from monthly 
temperature and salinity climatology and the thermal wind 
equations integrated upward from a predetermined "level 
of no motion”. The geostrophic currents are nondivergent . 
The sole contributor to the vertical advection component 
is, therefore, the instantaneous wind drift via the Ekman 
suction mechanism. 
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